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EXECUTIVE  SUMMARY 


A  polarization  modulated  analog  optical  link  is  presented  with  theoretical  treatment  of  nonlinearities  due 
to  polarization  distortion  loss  (PDL).  A  method  for  measuring  PDL  within  the  system  is  detailed.  A 
number  of  experiments  are  performed  which  characterize  the  modulator  nonlinearity  as  well  as  the 
interaction  of  multiple  nonlinearities.  The  primary  findings  include: 

•  Misalignment  of  the  modulation  arc  induced  second  order  nonlinearities  and  resembles  the 
negatives  effects  imposed  by  chromatic  dispersion. 

•  Partially  compensating  for  PDL  distortion  resulted  in  non-conformal  behavior  for  high  RF  power, 
indicating  that  multiple  nonlinearities  induce  second  order  distortion. 

•  Polarization  modulator  nonlinearity  contributes  to  second  order  distortion  and  indicated  that 
suppression  techniques  cannot  fully  mitigate  distortion  at  high  RF  power. 


E-1 


INVESTIGATIONS  OF  POLARIZATION  DEPENDENT  LOSS  IN  POLARIZATION 
MODULATED  ANALOG  OPTICAL  LINKS 


I  INTRODUCTION 

High  performance  analog  optical  links  for  military  applications  primarily  rely  on  the  most  common 
link  architecture  employing  intensity  modulation  with  direct  detection  (IMDD)  due  to  the  ease  of 
implementation.  Recently  a  special  type  of  phase  modulator  employing  TE  and  TM  modulation,  a 
polarization  modulator  (PolM),  has  gained  interest  in  the  field  of  microwave  photonics.  A  PolM  provides 
a  constant  intensity  format  that  when  aligned  properly  to  a  discriminator  can  function  analogous  to  a 
Mach-Zehnder  modulator  (MZM)  [1].  Many  works  have  detailed  the  use  of  a  PolM  followed  by  a 
polarization  controller  (PC)  which  can  be  used  to  align  to  a  polarization  beam  splitter  (PBS)  such  that  a 
photodetector  can  directly  demodulate  the  signal  [2-5],  but  have  not  detailed  an  algorithmic  way  in  which 
to  properly  align  the  modulation  arc  and  center  the  state  of  polarization  (SOP)  for  maximum  efficiency.  In 
links  employing  polarization  modulation,  various  processes  can  negatively  impact  performance  such  a 
polarization  dependent  loss  (PDL),  polarization  mode  dispersion  (PMD),  and  chromatic  dispersion  (CD). 
For  the  purpose  of  this  report  we  focus  on  the  impact  of  PDL  in  such  a  system.  Previously  we  have  shown 
through  first  aligning  a  PolM  for  intensity  demodulation  followed  by  the  addition  of  a  known  discrete 
source  of  PDL,  the  bias  point  of  the  link  can  be  adjusted  to  cancel  the  second  order  distortion  produced  by 
PDL  [6,  7]. 

This  report  will  expand  upon  the  work  in  [6]  and  [7],  where  first  a  general  description  of  PDL  in  the 
system  is  discussed  as  well  as  the  details  by  which  we  can  measure  and  calculate  the  PDL  through  system 
diagnostics.  Next  a  set  of  experiments  are  presented.  First  the  system  is  analyzed  for  the  ideal  and  a  non¬ 
ideal  alignment  when  no  PDF  is  present.  Next,  the  effect  of  including  a  discrete  PDF  and  only  partially 
compensating  the  resulting  phase  shift  is  presented.  Additionally,  when  no  PDF  is  present,  we  look  at  the 
difference  in  performance  when  adding  a  DC  bias  to  the  PolM  and  attempting  to  compensate  the  induced 
phase  shift  through  rotating  the  half-wave  plate  at  the  discriminator.  Finally  we  summarize  the  findings. 


2  POLARIZATION  DEPENDENT  LOSS:  THEORY  AND  SETUP 
2,1  GENERAL  THEORY  OF  PDL 


The  theoretical  description  relies  on  a  simple  model  of  PDL  as  an  ideal  partial  polarizer,  in  Jones 
space  [8]: 
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where  we  have  expanded  P  in  terms  of  the  2x2  identity  matrix,  (Tg  ,  and  the  Pauli  matrix  (Tj  .  For 
Hermitian  matrices,  as  expansion  in  the  Pauli  matrix  basix  set  results  in  4  numbers,  of  which  three  can  be 

associated  with  the  components  of  a  3-D  vector,  illustratively  P  ,  in  Stokes  space  [9]  -  [12].  In  the  simple 
particular  case  of  (1),  we  can  associate  the  operator 
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e,-{a,e,+a^e^+a^e^) 


=  p^cr^+PP-a  =  p^a^  +  P-a, 

with  average  field  transmission  /iQ=(l  +  t)/2  ,  pass  eigenstate  P  =  e^,  and  PDL  vector  P  =  PP , 
having  differential  field  magnitude  P  =  (l-t)/2  .We  note  for  a  general  polarizer  with  0  <  ^  <  1 , 


pl+p.pJ-±L 


pl-P-P  =  t 


(3) 


As  shown  earlier  [11],  [12],  once  the  operator  P  is  cast  in  form  (2),  it  can  be  generalized  to  any 

eigenstate  P .  Thus,  the  second  line  of  (2)  is  general  and  describes  any  orientation  for  the  pass  state  of  the 
link  PDL  or  the  receiver’s  analyzer. 

To  motivate  the  impact  of  PDL  impairments,  we  first  examine  the  transformation  of  state  |  incident 

on  a  partial  polarizer  to  its  output  state  |  ^ ,  in  terms  of  the  corresponding  Stokes  vectors  S  and  s ^  . 

Since 

=  (x|P'(7P|5)  (4) 

And  the  norm  of  ?  is  |5|  =  |  x)  =  ,  we  find 


Sf  =s 


\pI-P^)s  +  2{p„+P-s)P 


(5) 


Where  P^  =  P  P  and  S  is  the  unit  Stokes  vector  for  s  .  From  this,  the  power  transmission  is  found  as 

T  =  [pI+P^)  +  ^Po[P-s)  ,  (6) 

Giving  the  expected  transmission  characteristic  [16]. 

To  illustrate  the  impact  of  PDL  impairments,  we  examine  the  effect  of  P  on  states  of  polarization 
(SOPs)  for  several  value  of  pass  axis  transmission,?.  In  Fig.  1,  the  outer  spheres  represent  the 
unattenuated  SOPs  that  could  be  incident  on  the  partial  polarizer.  The  partial  polarizer  is  represented  by 
an  arbitrary  P :  In  this  case,  upward  to  the  left,  and  “behind”  the  page.  An  arc  of  modulated  SOPs  is 
shown  on  the  “far  side”  of  the  sphere,  centered  on  the  quiescent  SOP,  Xg  at  90°  to  the  x  (or  gj  )  vector  in 
this  Stokes/Poincare  representation,  as  in  Fig.  1(b).  Such  an  arc  of  modulated  SOPs  would  therefore  be 
ideally  biased.  PDL,  characterized  by  P ,  distorts  this  ideal  operation  as  indicated  by  the  inner  surfaces  in 
each  of  the  sub-figures,  drawn  to  represent  the  transmitted  power,  T  ,  for  each  SOP,  as  in  (6).  Four  such 
surfaces,  corresponding  to  different  values  of  t  in  (1)  are  shown,  with  (a)-(d)  running  from  0.9  to  0.1. 
Also  shown  are  the  arcs,  inside  those  surfaces,  corresponding  to  SOPs  emerging  from  the  partial 
polarizer.  For  t  =  0.9,  Fig.  1(a),  the  input  and  output  arcs  are  similar.  For  t  =  0.7,  Fig.  1(c),  the  presence 
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of  a  “dimple”  in  the  rejected  SOP  (opposite  to  P  )  is  becoming  evident  [8],  as  well  as  indications  that  the 
arc  is  distorted  in  direction.  For  high  PDL  t  =  \.  Fig.  1(d),  a  dimple  representing  the  rejected  SOP’s 
transmission  is  visible,  and  the  transmitted  SOPs  are  are  nearly  coincident  with  pass  axis  line. 

We  model  the  link  as  launch  state  evolving  under  the  operators  for  the  modulator  (M)  ,  PDL 

(P)  ,  and  polarization  discriminator  (D)  operators.  At  the  photodetector,  output  state  is  given  by 

|x^.)  =  DPM|xo)  (7) 

And  the  normalized  optical  intensity  is  found,  through  squaring  its  norm,  as 

~  {s  f  I  ‘^/  )  =  {^0 1  M^P^D^D  P  M  I  Xq  ) .  (8) 


P  and  D  are  covered  by  (l)-(3),  while  M  is  modeled  by 
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Which  introduces  a  phase  difference  of  (f)  between  the  x  and  y  components  of  the  field.  This  can  be 
expressed  as 


M„  ((!>)  = 


cos 


[(j) !  1)-  j  sm[(l) !  2) 
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0  cos(^zi/2)  +  7sin(^2i/2) 

=  cos(^zi/2)a'o  -  jsm{^ij)l  2)M -a 


(10) 


where  M  =  represents  the  slow  eigenstate  of  the  modulator.  Again,  a  similarity  transform  allows  this 
to  be  generalized  to  any  M  [1 1],[12]. 

Using  form  (3),  and  employing  (10)  for  operator  M ,  we  substitute  (8)  to  solve  for  the  photocurrent. 
Assuming  a  negligible  dispersion  and  an  ideal  modulator,  the  detector  photocurrent  is  given  by 


/jgt  (^z^)  =  — (ci +  AxM  ■  s -{^A- sin^zi 


cos^zi 


(11) 


Where  we  have  defined 


A  =  CjP  +  c^D 


(12) 


with 


Cl  =[dl+D-D')i^pl  +P-P)  +  4dQPoP-D 
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(13) 


Cj  =  2/>q  ^  Jq  +  D  •  D  j  +  AcIqP  ■  D 

c^=2d^[pl-P-P)  . 

Eqn.  (11)  is  the  central  theoretical  result  of  our  investigation.  It  is  a  description  of  the  output 

photocurrent  for  a  link  with  discrete  PDL  characterized  by  P  (but  otherwise  ideal),  expressed  as  a 
function  of  the  applied  modulator  phase,  (f)  ,  and  the  link  parameters. 

Before  proceeding  to  the  link  analysis,  we  make  contact  with  earlier  work.  For  the  special  case  of  no 
modulator,  ^  =  0  and  M  =  0 ,  and  (11)  reduces  to 

W=Ci  +  A-s.  (14) 

This  photocurrent  is  equivalent  to  the  power  transmission  T  and  agrees  with  the  general  form  in  Ref 
[13].  Furthermore,  Cj  is  an  explicit  form  for  the  reference’s  Fqn.  (5),  and  A  is  an  explicit  form  for  that 

reference’s  Fqn.  (6).  Thus,  this  special  case  of  (11)  agrees  with  the  earlier  transmission  form  and 
concatenation  rule  [13]. 

Returning  to  the  link  equations,  (11),  what  is  perhaps  most  surprising  is  that  it  has  the  general  form 

ijet  =  Cj  +  Cj  sin  ^  +  C3  COS  ^  (15) 


With  the  constants  C.  are  defined  explicitly  in  the  equation.  The  last  two  terms  are  somewhat  surprising. 
They  imply  that  the  time-varying  signal,  the  part  dependent  on  the  modulator  phase  (j)  ,  is  quite  similar  to 
the  ideal  behavior  in  a  polarization  modulated  link.  That  is,  the  output  in  (11)  is  a  sinusoidal  variation  in 

(j)  with  amplitude  and  a  phase  offset  of  tan  ‘  (C3  /  C2)  .  In  other  words,  apart  from  the  loss 

of  signal,  after  accounting  for  all  the  distortions  due  to  unintentional  amplitude  modulation  and  to  arc  re¬ 
orientation  in  Fig.  1,  the  output  has  exactly  the  same  form  as  an  ideal  interferometer  which  happens  to  be 
biased  off  quadrature  by  tan  '  (C3  /Cj)  .  To  put  this  in  context  of  Fig.  1,  the  projection  onto  the  X  for 

SOPs  along  each  of  the  four  arcs  (i.e.,  using  great  circle  arcs  connecting  each  point  to  .X )  is  the  same  as 
the  projection  of  the  SOPs  along  the  initial  arc,  after  a  slight  displacement  along  the  equator  and  an 
overall  attenuation.  We  will  experimentally  demonstrate  this  counter-intuitive  result. 

Since  M  -  M  =  \  for  our  assumed  ideal  modulator,  we  can  simplify  (11): 


W  =Ci+l 


{A- M^i^M  -8^+  AxM -s -{^A- M^i^M -s^  sin^  +  2-?cos^zi.  (16) 


Further  assuming  a  perfect  polarizer  at  the  receiver,  it  can  be  shown  that  (13)  reduces  to 
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Which  leads  to  a  direct  expression  for  (12): 


A 


o  .  2  O 

cos  —  +  /sin  — 

[2  2  J 

1  2  J 

(18) 


Where  O  is  the  angle  measure  of  the  great 


2.2  POLARIZATION  MODULATION  LINK  SETUP 


Fig.  2.1:  Experimental  Setup.  Laser  is  polarization  modulated  (PolM),  andlaunched  into  a  short  optical  link.  In  ideal 
link,  the  modulated  signal  is  developed  with  polarization  beam  splitter  (PBS),  detected  with  photodiode  (PD),  and 
the  RF  harmonics  are  analyzed  with  electrical  spectrum  analyzer  (ESA).  Polarization-dependent  loss  is  introduced  in 
a  test  jig  (PDL).  The  relationship  between  the  SOPs  of  the  signal,  the  PBS,  and  the  PDL  are  analyzed  by  introducing 
well-defined  polarization  states  with  a  polarization  state  generator  (PSG)  and  analyzed  with  reference  to  a 
polarimeter  (Polar  Analyzer).  Rejected  optical  power  is  monitored  (Opt  PM)  for  diagnostics. 

We  constructed  a  basic  link  using  a  Versawave  40  Gb/s  polarization  modulator  [14]  that  permitted 
both  RF  modulation  and  DC  bias  control.  The  light  from  a  DFB  laser,  after  modulation,  passed  through  a 
coupler,  was  amplified,  passed  through  a  set  of  polarization  controllers  (PCI A,  PCI)  to  a  tap  coupler. 
(We  use  two  polarization  controllers  for  experimental  convenience:  in  effect  we  must  set  not  only  ~s0  on 
the  correct  point  of  the  equator,  but  must  also  put  the  entire  modulated  arc  on  the  equator  [15].)  At  that 
tap  coupler,  a  portion  of  the  light  was  sent  to  a  polarimeter  (Polar  Analyzer)  which  we  used  to  monitor  the 
experiment,  as  detailed  below.  Light  from  the  other  port  was  sent  through  another  polarization  controller 
(PC2)  and  then  into  a  bulk  beam  expander  and  refocused  back  to  fiber.  The  collimated  beam  in  this  beam 
expander  was  used  to  introduce  the  various  PDL  elements,  as  described  below.  For  convenience,  a  half¬ 
wave  plate  (HWP)  on  a  rotatable  stage  was  mounted  in  a  second  bulk  optic  beam  expander.  This  provided 
a  reproducible  way  to  apply  and  invert  a  set  of  polarization  transformations,  also  described  below.  In  ideal 
operation,  the  net  result  of  the  polarization  evolution  from  the  PolM  output,  through  the  link  and  HWP,  to 
a  polarization  analyzer  (bulk  polarization  beam  splitting  cube  (PBS)),  was  that  the  light  at  the  PBS  input 
was  polarized  at  quadrature,  and  the  detected  light  at  the  photodetector  (PD)  mimicked  an  ideal 
interferometric  link  [16].  For  diagnostic  purposes,  the  power  at  the  rejected  PBS  port  and  the  DC 
photocurrent  on  the  output  photodiode  were  recorded.  An  electrical  spectrum  analyzer  (ESA)  was  used  to 
monitor  the  first  4  harmonics  of  the  detected  power  spectrum  as  a  means  of  quantifying  PDL  impairments 
introduced  into  the  link. 
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To  set  up  the  idealized  link,  we  followed  a  proeedure  established  earlier  [15]  that  ensures  ideal 
operation  by  reprodueibly  orienting  the  modulated  SOPs  which  enter  the  PBS.  The  link  is  then  biased  at 
quadrature,  there  is  no  PDL  in  the  test  jig,  and  the  SOPs  corresponding  to  the  modulator  eigenstate  ~M 
and  the  PBS  pass  axis,  ~D  have  images  which  are  marked  on  the  polarimeter  in  Fig.  2.3. 

2.3  PDL  MEASUREMENT  ROUTINE 

Once  the  operation  of  the  link  had  been  established,  we  then  introduced  PDL  by  dropping  in  a  home¬ 
made  test  jig  which  held  two  uncoated  optical  flats  held  at  various  fixed  angles.  An  angle-dependent  PDL 
is  introduced  by  tilting  since  the  s-  and  p-polarized  components  striking  both  surfaces  of  a  flat  experience 
different  Fresnel  reflection  coefficients  [8].  We  used  two  identical  fiats  oriented  with  opposite  angles, 
similar  to  an  earlier  description  [17],  so  that  the  PDL  is  doubled  while  walkoff  through  the  beam  expander 
is  cancelled.  To  measure  the  added  PDL,  the  main  laser  is  shut  off  A  polarization  state  generator  (PSG) 
from  General  Photonics  is  used  to  generate  a  test  signal  that  can  scan  through  the  6  degenerate  states.  At 
each  generated  state  the  SOP  at  the  PA  is  recorded  as  well  as  the  optical  power  received  at  the  power 
meter  at  the  1%  tap  of  the  coupler  placed  before  the  FIWP.  The  measurement  routine  allows  for  the 
determination  of  the  pass  axis  represented  at  the  PA  as  well  as  the  total  PDL  measured  in  the  link  through 
a  Mathematica  routine  employing  a  least  squares  fit  (Appendix  A).  The  measured  PDL  can  then  be 
compared  against  the  expected  PDL  determined  by  the  theoretical  calculations  based  on  the  angle  of 
rotation  of  the  optical  fiats.  Ideally  a  system  with  no  PDL  will  record  constant  optical  power  when 
performing  a  sweep  with  the  PSG.  As  illustrated  in  Fig.  2.3,  a  shaded  surface  is  shown  inside  the  Poincare 
sphere  which  represents  a  discrete  amount  of  PDL,  where  a  dimple  can  be  seen  that  is  orthogonal  to  the 
pass  axis.  Thus  one  can  conceptualize  how  the  measurement  of  6  points  on  the  shaded  surface  can  be 
translated  into  the  approximate  value  of  the  total  PDL  as  well  as  delineate  the  pass  axis,  all  from  the 
simple  routine.  Conceivably  a  well-designed  system  could  incorporate  a  built-in  PDL  measurement 
function  to  measure  the  total  PDL.  It  should  be  noted  that  this  routine  has  only  been  used  for  one  discrete 
source  of  PDL,  when  in  fact  the  total  system  may  contain  more  than  one  source  of  PDL.  The  theoretical 
development  to  handle  such  a  scenario  is  ongoing. 


Fig.  2.3:  Shaded  surface  graphs  60°  PDL  with  a  pass  axis  misaligned  with  the  modulator  pass  axis.  Orthogonal  to 
the  pass  axis  a  dimple  is  observed  in  the  PDL  surface.  The  surface  represents  the  non-uniform  optical  power  loss  at 
all  SOPs.  The  black  arc  represents  the  ideal  modulation  arc  along  the  great  surface  when  no  PDL  is  present.  The 
altered  modulation  arc  is  shown  in  magenta  where  the  arc  is  inside  the  Poincare  sphere  and  shifted  from  the  black 
arc. 
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Once  a  given  PDL  was  introduced  and  the  diagnostic  routine  described  above  was  performed,  we 
took  a  series  of  “characterization  curves”  to  quantify  the  effect  of  the  impairment.  The  characterization 
curves  for  each  of  the  PDL  levels  were  taken  by  applying  sinusoidal  2  GHz  signals  (at  a  sequence  of  RF 
powers)  to  the  polarization  modulator’s  input  port.  Each  RF  power  corresponds  to  a  well-defined  voltage 
excursion  amplitude,  which  in  turn  corresponds  to  a  well-defined  arc  length  of  SOPs.  For  each  such  input 
RF  power  level,  we  recorded  and  plotted  the  first  four  harmonic  RF  powers  from  the  ESA.  Even  for  ideal 
operation,  one  expects  the  transfer  function  (1)  to  generate  odd  harmonics  of  the  2  GHz  fundamental, 
while  that  quadrature  bias  suppresses  all  even  harmonics.  We  found  that  to  fully  characterize  the  outputs, 
all  4  harmonics  were  needed.  In  a  separate  experiment,  the  frequency-dependent  RF  losses  of  our  coax 
cables  were  measured  at  all  4  harmonics  to  correct  the  raw  data  from  the  ESA.  The  details  of 
measurements  for  PDE  with  optical  flats  at  various  angles  the  corresponding  theoretical  curves  can  be 
found  in  [6].  As  detailed  in  the  theoretical  treatment  of  PDE  in  [6],  PDE  corresponds  to  a  phase  shift  in 
the  link  which  manifests  as  even-order  distortion.  Consequently,  we  detailed  a  procedure  for  mitigating 
the  incurred  even-order  distortion  through  an  appropriately  applied  DC  bias  voltage  at  the  PolM  or 
equivalently  by  rotating  the  HWP  [7].  The  primary  focus  of  this  report  is  to  detail  a  series  of  experiments 
performed  where  we  interrogated  various  performance  aspects  of  the  link  in  Fig.  2.1. 


3  MISALIGNING  THE  MODULATION  ARC 

Up  to  this  point  the  system  has  been  operated  by  first  performing  the  alignment  procedure,  where  the 
center  of  modulation  is  located  at  45  degrees  to  the  pass  axis  and  the  modulation  arc  resides  along  the 
great  circle  as  shown  in  solid  red  in  Fig.  3.1.  Previous  work  has  described  the  effects  of  chromatic 
dispersion  on  a  polarization  modulated  link,  where  the  nonlinearity  manifests  as  a  misalignment  in  the  arc 
[18].  We  can  simulate  a  similar  case  experimentally  by  purposefully  misaligning  the  arc  during  our 
procedure,  represented  as  shown  in  the  dotted  curve  in  Fig.  3.1.  As  with  the  polarization  alignment 
procedure  detailed  previously,  the  SOP  is  aligned  to  the  top  of  the  sphere  and  represents  the  middle  of  the 
modulation  arc  (seen  in  Fig.  3.1  as  the  black  dot)  [1].  The  arc  must  then  pass  through  what  is  referred  to 
as  point  B  which  will  lie  somewhere  on  the  great  circle  when  using  the  PA  as  the  visualizer.  The 
calibration  procedure  is  first  performed  with  proper  alignment  as  described.  The  resulting  fundamental 
and  2"'*  through  4*  harmonics  for  a  2  GHz  input  signal  as  a  function  of  input  RF  power  are  shown  in  black 
in  Fig.  3.2.  Next  the  alignment  is  adjusted  such  that  a  new  B  is  recorded,  designated  B’,  which  is 
misaligned  by  10°.  Again  the  fundamental,  2"'*  and  4*  harmonics  are  recorded  and  plotted  in  Fig.  3.2  in 
blue,  red  and  green  respectively.  Although  the  fundamental  and  4*  harmonic  remain  the  same,  the  2““* 
harmonic  increases  significantly  at  low  input  RF  power.  As  has  been  shown  previously,  chromatic 
dispersion  in  an  intensity  modulated  link  results  in  increased  second  order  distortion  [19].  One  could  infer 
from  this  experiment  that  a  well-controlled  rotation  of  the  modulation  arc  about  the  axis  would  allow  one 
to  probe  the  effects  of  chromatic  dispersion  in  polarization  modulated  links  more  thoroughly  with  a 
combination  of  misaligned  arc  measurements  and  appropriate  lengths  of  fiber  to  validate  the  results. 
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Fig.  3.1:  Solid  red  arc  shows  ideal  alignment  at  input  to  PBS.  Arc  is  purposefully  misaligned  (red  dotted) 
to  investigate  resulting  nonlinearities. 


Fig.  3.2:  Fundamental,  2"‘*  and  4*  harmonics  for  2  GHz  input  signal  for  ideal  aligned  case  (black)  and  with  arc 
misaligned  by  10°  (color). 
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4  PARTIALLY  COMPENSATING  PDL  WITH  DC  BIAS 


As  was  theoretically  and  experimentally  proven,  PDL  generated  even-order  distortion  can  be  mitigated 
by  properly  adjusting  the  bias  phase  through  either  the  DC  bias  or  HWP  [7].  In  the  next  set  of 
experiments  we  look  at  the  effects  on  distortion  for  when  the  bias  phase  is  not  adjusted  to  the  proper 
cancellation  condition.  In  all  the  figures,  the  ideal  data  is  shown  in  black,  where  the  bias  voltage  was  4. 1 
V  for  a  PDL  of  ~  3dB.  The  colored  data  in  Fig.  4.1  shows  the  PDL  induced  even  order  distortion, 
previously  documented  with  the  theory  described  in  section  2,  evident  when  no  DC  bias  is  applied  and  the 
HWP  is  left  unaltered  [6].  Using  this  as  the  reference  point,  in  theory  the  curves  should  be  between  the 
red  and  black  points  for  the  2'“^  harmonic  shown  in  Fig.  4.1  and  the  green  and  black  points  for  the  4* 
harmonic.  In  Fig.  4.2,  3  V  is  applied  to  the  DC  bias  port,  and  as  expected  both  the  red  and  green  data 
move  closer  to  the  ideal  black  data.  Finally,  the  PDL  is  over  compensated  by  applying  5  V  (0.9  V  more 
than  the  ideal  condition  of  4.1V).  The  red  and  green  curves  for  low  RF  power  do  indeed  behave  as 
expected,  where  they  are  higher  than  the  ideal  data  in  black.  However  at  higher  RF  power  (  >  4dBm)  the 
red  and  green  data  are  below  the  black  data.  When  considering  only  one  source  of  nonlinearity,  the 
second  order  distortion  will  behave  as  a  slope  of  two.  However  when  the  model  incorporates  multiple 
nonlinearities  that  may  all  have  different  phases,  such  as  due  to  the  PolM,  the  photodiode  and  polarization 
dependent  loss,  the  second  order  distortion  at  the  output  becomes  much  more  complicated.  The  model  we 
have  built  in  Mathematica  (Appendix  B),  allows  us  to  get  a  better  picture  of  the  different  distortions  that 
all  interact  within  the  system.  However  in  order  to  tease  out  each  nonlinearity  a  set  of  experiments  needs 
to  be  determined  to  separate  nonlinearities  that  are  all  on  the  same  order  of  magnitude.  If  one  nonlinearity 
is  dominant,  the  model  can  easily  predict  the  curves  as  was  the  case  in  our  previous  work  [6].  The 
additional  experiments  to  more  thoroughly  investigate  the  model  are  outside  the  scope  of  this  report  and 
will  be  left  for  future  research. 


Fig.  4.1:  PDL  of  ~  3  dB  (optical  flats  at  60  degree)  is  added  to  the  system.  Fundamental,  2"^*  and  4*  harmonics  are 
plotted  for  2  GHz  input.  System  is  aligned  and  HWP  is  rotated  to  minimize  order  harmonic  at  -16  dBm  input  RF 
power,  after  which  the  PDL  is  added  (Color).  Same  procedure  as  for  data  shown  in  color,  but  an  additional  4.1  V 
bias  is  applied  to  the  DC  input  of  the  PolM  to  minimize  the  2"‘*  harmonic  at  -16  dBm  input  RF  power  after  adding 


PDL  (Black). 


9 


Fig.  4.2:  Curves  for  4.1  V  bias  applied  after  PDL  added  as  seen  in  Fib  6a  (Blaek).  Experiment  is  repeated  but  only 
3V  of  bias  is  applied  to  partially  eompensate  2"^*  harmonie  distortion  due  to  PDL  (Color). 


Input  RF  Power  (dBm) 

Fig.  4.3:  Curves  for  4.1  V  bias  applied  after  PDL  added  as  seen  in  Fib  6a  (Blaek).  Experiment  is  repeated  with  5  V 
bias  applied  to  over  eompensate  2nd  harmonie  distortion  due  to  PDL  (Color). 
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5  INVESTIGATING  NONLINEARITY  IN  A  POLARIZATION  MODULATOR 


Lastly,  a  set  of  experiments  were  eonducted  to  investigate  the  nonlinearity  due  to  the  PolM.  It  has 
been  assumed  up  to  this  point  that  the  DC  bias  at  the  PolM  and  the  rotating  HWP  are  equivalent  in  terms 
of  bias  phase  eontrol  points.  As  was  shown  previously,  the  DC  bias  exhibits  its  own  nonlinear  optieal 
transfer  funetion  that  may  impose  sueh  nonlinearity  on  the  system.  The  link  was  setup  and  aligned  as 
before  initially  with  no  DC  bias  applied.  No  PDL  is  added  and  the  2"'^  harmonic  distortion  is  mitigated  by 
adjusting  the  HWP  properly.  In  the  following  figures,  the  resulting  data  from  the  fundamental,  2"“*  and  4* 
harmonics  is  plotted  in  black.  Next  we  applied  varying  DC  bias  voltages  to  the  PolM,  followed  by  the 
alignment  procedure.  Finally  the  HWP  is  again  rotated  in  order  to  mitigate  the  2"‘*  harmonic  distortion.  If 
adding  a  DC  bias  is  equivalent  to  adjusting  the  HWP,  the  output  of  both  experiments  should  be  the  same. 
Bias  voltages  of  3,  6,  8  and  10  V  were  tested  with  the  output  shown  in  Figs.  5. 1-5.4  respectively.  In  the 
first  case  3  V,  the  4*  harmonic  is  slightly  reduced  throughout,  while  the  2’“^  harmonic  is  the  same  at  low 
and  high  RF  power,  but  lower  in  the  mid  RF  power  range.  Thus  the  distortion  curves,  particularly  for  the 
2”“*  order  distortion  cannot  be  easily  described  by  a  slope  of  2,  as  is  generally  thought  of  for  2”'^  order 
distortion.  Clearly  other  nonlinearities  that  occur  within  the  system  dominate  and  potentially  more  than 
one,  as  the  photodiode  distortion  is  mostly  mitigated.  We  find  for  the  6  V  case  that  the  curve  shows  more 
distortion  at  the  low  RF  power  and  less  at  high  RF  power  for  the  2"'*  harmonic,  so  there  isn’t  a  clear 
benefit  to  using  either  the  HWP  or  the  DC  bias.  In  the  final  two  cases  for  8  V  and  10  V  the  2"‘^  harmonic 
again  is  similar  at  lower  RF  power,  but  significantly  lower  at  high  RF  power,  upwards  of  1 5  dB  near  the  2 
dB  RF  input  point.  The  result  of  these  measurements  allows  us  to  conclude  that  the  distortion  generated  in 
the  link  does  not  behave  as  ideally  as  our  model  would  suggest,  and  that  the  bias  points  can  be  potentially 
used  in  tandem  to  mitigate  nonlinearity  even  further. 


Fig.  5.1:  Fundamental,  2nd  and  4th  harmonics  are  plotted  for  2  GHz  input.  System  is  aligned  and  HWP  is  rotated  to 
minimize  2nd  order  hannonic  at  -16  dBm  input  RF  power  (Black).  Same  procedure  as  for  data  shown  in  black 
except  before  alignment  is  performed  a  DC  bias  of  3  V  is  added  to  the  PolM  (Color). 
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Input  RF  Power  (dBm) 

Fig.  5.2:  Fundamental,  2nd  and  4th  harmonics  are  plotted  for  2  GFIz  input.  System  is  aligned  and  FIWP  is 
rotated  to  minimize  2nd  order  harmonic  at  -16  dBm  input  RF  power  (Black).  Same  procedure  as  for  data 
shown  in  black  except  before  alignment  is  performed  a  DC  bias  of  6  V  is  added  to  the  PolM  (Color). 


Fig.  5.3:  Fundamental,  2nd  and  4th  harmonics  are  plotted  for  2  GFIz  input.  System  is  aligned  and  FIWP  is 
rotated  to  minimize  2nd  order  harmonic  at  -16  dBm  input  RF  power  (Black).  Same  procedure  as  for  data 
shown  in  black  except  before  alignment  is  performed  a  DC  bias  of  8  V  is  added  to  the  PolM  (Color). 
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Fig.  5.4:  Fundamental,  2nd  and  4th  harmonics  are  plotted  for  2  GFIz  input.  System  is  aligned  and  FIWP  is 
rotated  to  minimize  2nd  order  harmonic  at  -16  dBm  input  RF  power  (Black).  Same  procedure  as  for  data 
shown  in  black  except  before  alignment  is  performed  a  DC  bias  of  lOV  is  added  to  the  PolM  (Color). 
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6  SUMMARY  AND  CONCLUSIONS 

In  summary  the  theory  for  a  discrete  PDF  in  a  polarization  modulated  link  was  detailed.  The  link  was 
presented  along  with  the  measurement  routine  needed  to  measure  PDF  in  the  system.  The  link  alignment 
for  optimal  intensity  modulation  at  the  output  is  done  before  PDF  is  added  to  the  system,  which  does  not 
represent  a  real  world  case  where  unknown  PDF  would  be  part  of  the  system  before  alignment.  None  the 
less  the  link  was  characterized  for  cases  when  the  modulation  arc  is  slightly  misaligned,  which  shows 
contributed  second  order  nonlinearity  and  the  experiment  is  similar  to  the  arc  rotation  that  occurs  with 
chromatic  dispersion.  A  set  of  experiments  were  performed  to  investigate  partially  compensating  the 
PDF.  The  results  indicate  that  multiple  nonlinearities  are  acting  within  the  link  and  their  additive 
magnitudes  and  phases  provide  a  much  more  complex  result  than  a  system  with  a  single  nonlinearity  that 
dominates.  Thus  for  mitigation  purposes,  we  must  be  careful  to  account  for  all  sources  of  nonlinearity  as 
suppressing  the  largest  will  leave  other  smaller  nonlinearities  in  the  system.  Finally,  the  nonlinearity  of 
the  PolM  was  investigated.  We  find  that  when  applying  a  DC  bias,  the  distortion  at  higher  power  behaves 
much  differently  and  cannot  be  completely  mitigated  with  the  FIWP.  This  echoes  the  conclusion  from  the 
partial  PDF  compensation  experiments,  that  at  higher  RF  power,  with  multiple  nonlinearities  in  the 
system,  suppression  techniques  cannot  fully  mitigate  strong  distortion. 
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APPENDIX  A 


PDL  Inverse  Routine 


dB2lin[<lB_]  :=  10*  (dB/ 10);  (*  convert  dB  to  linear  units  *) 

unit[s_]  :=  If[Nonn[s]  =0,  (0,  0,  1),  s/Nonn[s]  ]  (*  make  3D  unit  vector 

(*  power  of  ’sop’  after  going  thru  polarizer  oriented  at  'pop*  where 
pi  is  FIELD  xmsn  thru  'pop'  and  p2  is  FIELD  xmsn  thru  '-pop'  *) 
pwr  [sop_,  pop_,  pl_,  p2_J  :=  (»  power  *) 

((pl*2+p2*2)  +  (pl*2  -  p2*2)  (unit ( sop]  .unit [pop] ))/ 2; 

f indpdlLIN  [polpwrs  ]  :=  Module  [{  sdat ,  soln,  ppop,  pdlfit,  tran. 


yy,  zz,  xp,  yp,  zp,  pfx,  pfy}, 

( *  OSAGE  •  ) 

(*  polpwrs  is  set  of  {SOPi ,  P(lin) [points  for  6 
points  ...  *} 

(«  {  {3ll,sl2,sl3,pl) , (s21,s22,s23,p2} ,  *> 

(•  ...  { s61 , 362 , s63 ,p6}  }  *) 

<*  subroutines  needed:  pwr,  \init,dB21in  *} 

(*  FindPit  [data,  f  [pars ,'vars ]  ,  pars,  vars]  *) 


soln  =  FindFit  [polpwrs,  pwr[{xx,  yy,  zz>,  {xp,  yp,  zp} ,  pfx,  pfy], 
(xp,  yp,  zp,  pfx,  pfy},  {XX,  yy,  zz},  AccuracyGoal -►  4]  ; 
ppop  =  unit[{xp,  yp,  zp}  /.  soln];  <*  estimate  for  pass  axis  •) 

(•  once  FIELD  xmsns  are  found,  square  for  intensity  *) 
pdlfit  =  10  Logic  [  (pfx  /  pfy)  *2]  /.  soln; 

(•  estimate  for  PDL  magnitude  •) 

(pdlfit,  ppop}  («  return  fit  *) 

]; 


dataSOP  =  { 


{+  .7668, 

+ .5612, 

- .3115, 

{+ .5801, 

- .8135, 

- .0257, 

{- .5615, 

+ .8235, 

+ .0763, 

1 

GO 

00 

- .5531, 

+ .3654, 

{+  .3055, 

+ .1574, 

+ .9394, 

{  - .2862, 

- .1773, 

- .9415, 

}; 

f indpdlLIN [dataSOP] 

{3.1654,  {-0.13'19e,  0.747118 


57.12},  (•  SOP  powers  in  uW  •) 

41.31}, 

69.07} , 

53.52}  , 

68.58}  , 

41.77} 


0.650842} } 
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APPENDIX  B 


Lab  data  fitting  140801  60v2 


Functions  and  constants 


ln[1]:=  dBtop  =  0; 

dBbot  =  -140; 

dBlft  =  -25; 

dBrt  =  15; 

pltmg  =  ({dBlft,  dBrt},  (dBbot,  dBtop}}; 
color  =  (Blue,  Red,  Green,  Brown}; 
npts  =  100 ; 

usin  =  Table[Sin[2  Trn  .  npts] ,  (n,  0,npts-l}l; 

; •  unit  sine  function  for  sod  input  • ) 

corr2  =  1.15;  (♦  TBESK  ARE  COAX  LOSS  CORHECTICWS  FOR  2, 4, 6, 6  GHz 

corr4  =  1.62; 

corr6  =  1.92; 

corrO  =  2.46; 
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In[13]:= 


a2dB[a_]  :=  10  Logl0[0.5  a'^2] ;  {•  convert  sine  aiipl  to  power  «) 

dB2a[d_]  :=  Sqrt[2  x  lO*  (d  /  10)] ;  («  convert  power  to  anplitude  •) 

3pec[power_]  :=  Module  [{ft,  spectrxna},  f*  calc  1st  4  Hannonlc  powers  ») 
ft  =  Fourier [power] ; 
spectruK  =  4  »  ( .  5  Abs  [ft]  *2)  /npts; 

{spectruii[  [2]  ] ,  spectrum]  [3]  ] ,  spectrum]  [4]  ] ,  spectrum]  [5]  ]  } 

]; 


{•  power  of  'sop'  after  going  thru  polarizer  oriented  at  'pop'  where 

pi  is  FIELD  zmsn  thru  'pop'  and  p2  is  FIELD  xmsn  thru  '-pop'  •) 

pwr[sc>p_,  ppp_,  pl_,  p2_]  :=  ((pl''2  +  p2''2)  *  (pl*2  -  p2'^2)  <unit[sop]  .unit[pop]))  '  2 
unit[s_]  :=  If[Norm[s]  =  0,  {0,  0,  1},  s  'Norm[s]];  (*  make  3D  unit  vector  •) 
dB21in[dB_]  :=  10'' (dB  /  10)  ;  (*  convert  dB  to  linear  units  *) 

lin2dB[lin_]  :=  10  LoglO  [Iln]  ;  (»  convert  linear  units  to  dB  *) 
findpdlLIN[poIi>wrs_]  :=  Module [{sdat,  soln,  ppop,  pdlfit,  tran, 
yy»  zz»  zpr  yp,  zp,  pfx,  pfy}, 

( •  OSAGE  . ) 

(•  polpwrs  is  set  of  {SQPj^,  P<lin)  }points  for  6  points  ...  •) 

(•  (  (sll,sl2,sl3,pl} ,  {s21,s22,s23,p2) ,  ...  (s61 ,s62 ,363,p6]  }  .) 

(•  subroutines  needed:  pwr,  unit,  dB21in  •) 

(•  FindFit[data,  f [pars;vars] ,  pars,  vars]  •) 

soln  =  FindFit [polpwrs,  pwr[{xx,  yy,  zz),  {xp,  yp,  zp},  pfx,  pfy] , 

{zp,  yp,  zp,  pfx,  pfy) ,  {XX,  yy,  zz}]  ; 
ppop  =  unit[{xp,  yp,  zp}  /.  soln];  (*  estimate  for  pass  axis  «) 

:  •  once  FIELD  xmsns  are  found,  square  for  intensity  *) 

pdlfit  =  10  LoglO  [  (pfx  /  pfy) ''2]  /.  soln;  (•  estimate  for  PDL  magnitude  •) 
(ixilfit,  ppop} 

]; 

correct [dat3nat_]  :=  Module[{tr},  (»  Apply  coax  loss  corrections  to  RF  data  ») 
tr  =  Transpose  [datut] ; 

Transpose [{tr[  [1]  ] ,  tr] [2] ]  +  corr2,  tr [ [3] ]  +  corr4,  tr [ [4] ]  ♦  corr6, 
tr[[5]]  tcorrO}] 

]; 

unit[s_]  :=  If[Norm[s]  =0,  {0,  0,  1},  s/Norm[s]];  (*  makes  \mit  vector*) 
flip[r_,  s_]  :=  2unit[r]  (unit[r].s)  -  s;  (•  to  flip  180,  s  spar  -  sperp  •) 
flipax  [dd_,  ss_]  :=  unit  [unit  [  cid]  +  unit]  ss]  ] ;  (•  parallelogram  •> 


18 


Data  entry  for  RF  mn  and  6  point  SOP  meas  for  PDL  vector 


ln[25];= 

datavO  =  {  (•THIS  IS  RAH  DATA  *) 


{-20, 

-48.12, 

-123.57, 

-135.37,  - 

-135.16,  2.106,  2.354}, 

{-18, 

-46.08, 

-119.819, 

-131.51, 

-137.43,  2.103,  2.348}, 

{-16, 

-44.09, 

-115.089, 

-127.00, 

-136.34,  2.103,  2.348}, 

{-14, 

-42.12, 

-111.311, 

-121.266 

,  -139.74,  2.109,  2.355} 

{-12, 

-40.15, 

-107.419, 

-115.79, 

-146.62,  2.106,  2.352}, 

{-10, 

-38.14, 

-104.49, 

-109.82,  • 

-145.81,  2.104,  2.346}, 

1 

o 

00 

-36.12, 

-98.989, 

-103.59, 

-141.21,  2.110,  2.345}, 

1 

o 

-34.13, 

-94.048, 

-97.62,  - 

135.76,  2.113,  2.339}, 

{-04, 

-32.14, 

-89.986, 

-91.55,  - 

151.93,  2.119,  2.335}, 

{-02, 

-30.19, 

-85.052, 

-85.59,  - 

141.1,  2.123,  2.324}, 

{+00, 

-28.18, 

-79.999, 

-79.54,  - 

127.38,  2.130,  2.314}, 

{  +  02, 

-26.29, 

-74.727, 

-73.66,  - 

123.78,  2.142,  2.300}, 

{  +  04, 

-24.39, 

-69.476, 

-67.63,  - 

116.337,  2.156,  2.269}, 

{+06, 

-22.50, 

-63.585, 

-61.74,  - 

107.711,  2.173,  2.228}, 

{  +  08, 

-20.77, 

-57.288, 

-55.81,  - 

98.059,  2.184,  2.176}, 

{  +  10, 

-19.15, 

-51.318, 

-50.01,  - 

88.410,  2.214,  2.096}, 

{  +  12, 

-17.84, 

-45.520, 

-44.30,  - 

78.295,  2.228,  1.989}}; 

dataSOP  =  ( 

{-.9795,  +.1822,  +.0908,  15.58},  {•  SOP  powers  in  uW  •) 

{+.3015,  +.8965,  -.3209,  12.83}, 

{-.1314,  -.9241,  +.3600,  13.12}, 

{+.9755,  -.2174,  -.0214,  10.70}, 

{+.2080,  +.3271,  +.9219,  16.60}, 

{-.0289,  -.3362,  -.9412,  9.343} 

}; 


datavl  =  correct [ da tavO ] ;  (•  ai^ly  RF  coax  losses  to  data  •) 
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Run  dau  through  the  model 


ln{291:=  sopout  =  findpdlLIH[dataSOP]  ; 

(•  returns  |extiiiction  in  cffi,  |pa5s  axis  p||  ») 

tp  =  10'*  (-sopout[  [1]  ]  ■  20)  ;  {*  txansnission  for  PDL  •) 

PPhatO  =  sopout [[2]]  ;  P  for  PDL  •) 

PPO  =  PPhatO  (l-tp>  ,  2  ;  <*  P  for  PDL  «) 
pO  =  (1  +  ip)  /2; 
pbar  =  (1  -  tp)  /  2 ; 

SS  =  {.9686,  -.2364,  .0610);  <•  from  experiaent  after  90  deg  rotation  •) 

DD  =  {0,  0,  1}; 

DDhat  =  unit[DD] ;  olaage  of  detector  axis  on  PA  *) 

MM  =  DDhat  «SS;  («  image  of  modulator  eigenstate  •) 

rr  =  flipax[DD,  SS] ;  find  HWP  flipping  axis  as  bisector  of  DD  and  SSrl 
PP  =  fliptrr,  PPO];  {•  flip  the  P  ■vector  image  on  the  PA  «) 

PPhat  =  unit[PP]; 

PHI  =  ArcCos [PPhat. unit[DD] ] ; 

AA  =  ((Cos[PHI  2])  *2  ♦  tp  {Sin[PHI  2])  *2)  PP  +  tp  DDhat  2; 

CC2  =-(MM«AA.SS+  (AA.MM)  (HM.SS)); 

CC3  =  AA.SS; 

ln[46];=  s3  =  S(irt[3]; 
thk  =  0.2; 

circle  =  Graphics  [{Thickness  [thX] ,  Circle [{0,  0)]}]; 
plus  = 

Graphics [{Thickness [thk],  Line[{{-.5,  0],  {0.5,  0),  {0,  0],  {0,  .5),  {0,  -.5)}]}]; 
triattgle  = 

Graphics [ 

{Thickness  [thk]  ,  Line  [{{-.5,  -s3  6],  {  .5 ,  -s3  6}  ,  {0 ,  s3  3}  ,  {- .5,  -s3  6))])]; 

square  = 

Graphics [ 

{Thickness [thk] ,  Line[{{-.5,  -.5),  {.5,  -.5),  {.5,  .5),  {-.5,  .5),  {-.5,  -.5})])]; 
tr  =  Transpose [da tavl] ;  t«  create  plots  of  data  points  to  help  fit  model  •) 
fundraw  =  Transpose[{tr [ [1] ] ,  tr[[2]])]; 
secondraw  = Transpose[{tr [ [1] ] ,  tr[[3]]}]; 
second trim  =  secondrav[ [3  ; ;  Length [ secondraw] ]] ; 
thirdraw  =  Transpose [ {tr[ [1] ] ,  tr[[4]]}]; 
thirdtrim  =  thirdraw[[3  ;;  Length [thirdraw] ]] ; 
fourthraw  = Transpose[ [tr [ [1] ] ,  tr[[5]]}]; 
fourthtrim  =  fourthraw[ [10  ;;  Length [fourthraw] ]] ; 

f\md  =  Lis tPlot[ fundraw,  PlotMarkers  -•  {circle,  .05},  PlotRange pltmg. 

Frame  -»  True,  Axes  False  ] ; 

third  =  Lis tPlot[ thirdtrim,  PlotMarkers  -•  {triangle,  .05}]  ; 
second  =  ListPlot [secondraw,  PlotMarkers  ->  (plus,  0.05}] ; 
second  =  ListPlot [secondtrim,  PlotMarkers  -•  (plus,  0.05}]; 
fourth  =  ListPlot  [fourthtrim,  PlotMarkers  [square,  .05}]; 


ln[65]:=  Transpose  [{ tr  [  [l]  ]  ,  tr[[2]]}]; 
ln[66]:=  Manipulate 
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Out(66]= 


K>ts  = 

Showe« 

(ListPlot[Transpose[pptsout]  [[»?]],  PlotStyle  {color[[fiF]]  ,  PointSize[.008]}, 
PlotRange  -» pltrng,  Fraoie  -»  True,  Axes  False  ]  &  @  {1,  2,  3,  4)) ; 

label  =  Graphics [{ Text [ Style ["PDL  60",  24,  FontFamily "Helvetica"],  {-20,  -20}]}] 
Sbow[ • •g3ext,g2ext,g4ext, • ) fund,  second,  third,  fourth,  ppts,  label, 

FraaeLabel  -*  ("PolM  input  dB*",  "BSA  output  dBm"}, 

LabelStyle  -.  Directive  [Large ,  FontFaiiily  "Helvetica"]]  , 

{{a2,  -.022},  -.4,  .4},  {{ib,  0.378},  -.6,  .6},  {{eps,  0.014},  -.4,  .4}, 

{{RFIH,  14.6},  10,  20},  {{RFODT,  2 
}  ,  0,  25}  ,  SaveDefinitions  -.  True 
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